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Based on the method proposed in [ H. S. Zong, W. M. Sun, Phys. Rev. D 78, 054001 (2008)], 
we calculate the equation of state (EOS) of QCD at zero temperature and finite quark chemical 
potential under the hard-dense-loop (HDL) approximation. A comparison between the EOS under 



,^ ■ HDL approximation and the cold, perturbative EOS of QCD proposed by Fraga, Pisarski and 

Q^. Schaffner-Bielich is made. It is found that the pressure under HDL approximation is generally 

^^ , smaller than the perturbative result. In addition, we also calculate the quark number suscepti 



bility (QNS) at finite temperature and finite chemical potential under hard-thermal/dense-loop 
(HTL/HDL) approximation and compare our results with the corresponding ones in the previous 
literature. 
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The investigation of the equation of state (EOS) for cold and dense strongly interacting 
matter and its consequence for the possible phases of quantum chromodynamics (QCD) 
plays a crucial role in the study of neutron stars in astrophysics [ll, l2| . At present lattice 
QCD provides the most powerful tool in studying QCD thermodynamics. While the Monte 
Carlo simulation works well and has achieved great success in dealing with finite temperature 
QCD thermodynamics, it encounters the notorious sign problem in the situation of finite 
quark chemical potential and cannot be applied [3l]. To circumvent the sign problem, several 
methods have been proposed. Among them are the reweighting technique ^, analytic con- 
;inuation from imaginary chemical potential |5|, |6| and Taylor expansion in fx around yU = 
7H9|. Although these methods have achieved some success in treating finite /i physics, they 
are far from being complete and even suffer from some limitations and drawbacks [10]. So 
we expect that continuum model studies should be complementary to the lattice simula- 
tions in the exploration of QCD thermodynamics at finite /i. At high enough temperature 
(T) and/or the quark chemical potential (/i), the hard-thermal/dense- loop (HTL/HDL) ap- 
proximation is thought to be a good approximation for QCD [HJ]. In the present paper, 
based on the method proposed in Ref. [12|, we try to calculate the EOS of QCD at zero 
temperature and high enough quark chemical potential under HDL approximation. Since 
the HDL approximation is a good approximation for QCD at high enough quark chemical 
potential, the EOS derived with the aid of it should have a relatively good QCD foundation 
and one expects that it can be applied to the study of neutron star. In addition, with the 
method proposed in the present paper we can calculate another important quantity for the 
study of QCD phase transition, the quark number susceptibility (QNS), in the HTL/HDL 
approximation. The QNS measures the intrinsic statistical fluctuations in a system close to 



thermal equilibrium and is thought to p 



point in the QCD phase diagram |13|-|20|. 



ay an important role in identifying the critical end 



Let us start with the renormalized partition function of QCD at zero T and finite /i which 
reads 

Z[jj] = / VquVquVAn exp \ -SR[qR, Qr, Ar] + / d'^x lJ,Z2qR{x)-f4qR{x) > , (1) 

where SR[qR,qR,AR\ is the standard renormalized Euclidean QCD action with qR being 
the renormalized quark fields with three flavors and three colors, Z2 = Z2{C^,A'^) is the 
quark wave-function renormalization constant {( is the renormalization point and A is the 



regularization mass-scale). Here we leave the ghost field term and its integration measure 
to be understood. 

The pressure density P(/i) is given by 

V{^^) = ^\nZ[^^l (2) 

where V is the four-volume normalizing factor. The above expression for 'P(/i) is just the 
EOS. From the pressure density one immediately obtains the quark number density: 

M = -^. (3) 



According to Ref. 12|], one has the following well-known result 

P(/^) = -f Tr{G«H74} = -iV.iV;^2/^tr{G^[/i](p)74}, (4) 

where G/j[/i](p) is the renormalized quark propagator at finite chemical potential and Nc, 
Nf denote the number of colors and fiavors, respectively. Here, Tr denotes trace operation 
over color, fiavor, Dirac and coordinate indices, while tr denotes trace operation over Dirac 
indices only. 

Integrating both sides of the equation p{fi) = dV{fi)/dfi gives 

^(/i) = n^^)\,=o + l'd^,'pi^^') = r{p\^,-N,NfZ^ j^ dp' j ^^tr {Gn[A{p)lA} , (5) 

where 'P{p)\ ^ is an integration constant which represents the pressure density at /i = 0. 
Here one should note that 'P{p)\ ^ is a constant independent of p and the whole nontrivial 
p dependence oiV{p) is contained in the integration term which is totally determined by the 
renormalized quark propagator Gfi[p]{p). Formally Eq. ([5]) provides a model- independent 
formula for calculating the pressure density at finite p. However, at present we still do not 
know how to calculate the quark propagator at finite chemical potential from first principles 
of QCD. So, when one uses Eq. ([5]) to calculate the pressure density, one has to resort to 
various QCD models or approximations. In this paper we will use the HDL a ppr oximation. 



The quark propagator under HTL/HDL approximation can be written as 11| 



^^^^) = -^^ Dj^) 2 ^^ DM 2 ' (^) 

where p = p/|p| and p4 = {2n -|- l)7rT in G Z) are the fermion Matsubara frequencies. The 
form of the functions D±{p) is 



2 

D±{p) = -ipA + p±\p\ + -^ 



tp4-p\ ( %Pi- P 



(7) 



where ruq = gJ(T'^ + jj? /'k'^)/Q is the quark thermal mass with g being the strong couphng 
constant, Qq and Qi ^-^e Legendre functions of the second kind. Substituting Eq. ([6]) into 
Eq. (j4]) gives 
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According to a famihar method in thermal field theory, the frequency sum in Eq. (J8]) can 
be calculated by means of contour integral 
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The integral contour is shown in Fig. [H Now let us calculate the integral. Defining 

/(P4) := 



1 



1 



- tanh — - 
2 2T 



and closing Ci and C2 with large half circles, one obtains the following 
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^^^^^^ = ^(-27r.)5:Res{/(C,)}, 



(9) 
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where (j are the poles of f{p4) located in the imaginary axis. The poles of 1/-D+ are 
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FIG. 1: The integral contour 



Ci = -i(w+ + /i) , Res< 
C,2 = i{uj- — /u) , Res ■ 






iZ. 



(12) 
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and the poles of 1/D_ are 



Ca = ?(w+ - /i) , Res 
C4 = —i{u- + /i) , Res • 
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where uj±{\'p\){> \p\) are the solutions of the equations 
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and Z-i-(IpI) are 
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From these results one can calculate the right hand side of Eq. (fTTI) : 
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The second term in the left hand side of Eq. (9) can be expressed as 



+ iZ^ tanh 



(19) 
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where p±(a;, \p\) are the familiar spectral functions of the quark propagator (a; = uj/\p\) 11| : 



p±(u;,ip1) = 2n[z^m)si^-^±m)) + z^mn^+^Ti\m] 
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Combining Eqs. (19), (20) and (9), one can obtain 
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where npiu)) is the Fermi distribution function 



1 



exp((X'/T) + 1 
Therefore the quark number density at finite /i and T is 

p(/i,r) 



(23) 
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From Eq. flM|) one can clearly see that the quark number density contains two terms, one 
is the contribution from the quasi-particle poles and the other is the contribution from the 
Landau damping. In the limit of T — )■ oo, p tends to the free quark gas result (Z+ — > 1, 
Z_ — 7- and the integral is dominated by p ~ a; ~ T region) . 

Before proceeding the calculation of EOS, let us turn to the calculation of the QNS, which 
is defined as 

dpip, T) 



x{f^,T) = 



dp, 



(25) 



As mentioned before, the QNS plays an important role in identifying the critical end point 
in the QCD phase diagram. There are many calculations of QNS under HTL approximation 
in the literature 2l|-|25|. Generally, in order to calculate the QNS, one has to calculate the 
quark propagator and the quark-meson vertex in the vector meson channel at zero total 
momentum separately. In Ref. ^], using the fact that d/dp can be replaced by d/d{—ip^) 
and with the help of the Ward-Takahashi identity, the authors calculate the QNS under 
HTL approximation only through the HTL quark propagator. Here a question arises: is the 
information of p dependence lost by such a method in Ref. J25|? Now we will answer this 
question in a straightforward way by calculating the QNS directly from Eqs. (!24|) and (!25ll . 



Because Z±, p± and u± depend on p only through niq, one would obtain the following 

dp{p,T,mg) dp{p,T,mg)dmg 
xif^,T) = ^- + TT- TT-- (26) 



dp 
Setting p = in the above equation gives 
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drrin 



dp 



x(/i = 0,T) 



dp 



+ 



fi=0 



dp{p, T, nig) drrig 
drrig dp 



(27) 
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One can easily find that the second term in Eq. (1271) is zero {{dmq/dij)^=Q 
QNS at finite T and zero yU is 



0) and the 



x(/i = 0,T) 



T 

2N,Nf 






Wj 



2r7i^ 



nF[uJ+)[l - nF{uj+)) + — 



-Q 



21711 



ni?((X'_)(l — npioj-)) 



T J (27r)3 
This is the same as the result in Ref. 



^ 77-1 [P+(w) + p^{uj)]nF{uj){l - np{uj)) 



M=0 



25| . as one expected in advance. Here it should 



be stressed that the method used to calculate the QNS in the present paper is difi^erent 
from the one used in Ref. (25| in which the Ward-Takahashi identity is adopted to avoid 
the differentiation over [i. The result obtained here can be regarded as a self-check of the 
calculation in Ref. [25|. As is pointed out in Ref. 25|], the QNS contains the contributions 
from both the Landau damping and the quasi-particle poles. This is different from the result 
in Ref. 23| (see Ref. 25| for more detail). 

After the discussion of QNS let us continue the calculation of EOS. The zero temperature 
result of p reads 



P(/^) 



2N,N^ 



(27r)3 






(29) 



{2TTy J-\p\ 2tt' 

It should be pointed out that according to Eq. (4), in order to calculate the quark number 
density one needs to know the quark propagator in the whole momentum range. However, 
the HTL/HDL approximation is only valid for external momentum much smaller than T or 
fi. Therefore, as in Ref. 25|, for a consistent calculation of the quark number density under 
the HDL approximation one should introduce a momentum cutoff A^dl below which the 
HDL quark propagator (6) is applicable. In our calculation we choose Ahdl = A'p/ti"- We 
make this choice because gfj^/n is an important energy scale to identify the soft momentum 
in HDL approximation (for HTL this would be gT, see Ref. (25|). If we limit the range of 
integration of the contribution of Landau damping in Eq. fl29|) to the region |p| < gfJ^/TT, we 
will get the numerical result of the quark number density shown in Fig. 2 {N^ = 3, Nj = 2). 
From Fig. 2 it can be seen that the quark number density under HDL approximation is 
slightly smaller than that of free quark gas (p/ree(p) = NcNffi^/Sn'^), which indicates that 
the attractive interaction dominate in the high density case. Here we want to stress that the 
result in Fig. 2 has minor sensitivity to the variation of the cut-off S'/x/tt. For example, if one 



(28) 



sets the cut-off to be 2gfi/'n', one finds that the change of p is about 10%. This is because 
the quarlc number density (and susceptibihty) is determined essentially by the quasi-particle 
contribution within the given approximations. 
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FIG. 2: The quark number density under HDL approximation 



Now we can calculate the EOS under HDL approximation (here we neglect the constant 
term V{fi)\^=Q because when applying our EOS to the study of neutron star, owing to the 
boundary condition imposed on the surface of neutron star, the constant term P(/i)|^=o does 
not contribute to the mass-radius relation) 



+N.N,[d,'j 



(2vr) 



[p4oj)+p4u)][l-29iuj-p')]. (30) 



{2tt)^ J-\p\ 27r 

The numerical results are shown in Fig. 3 with Pfree = NcNfp"^/ {1271"^) being the pressure 
density of the free quark gas. It can be seen that the pressure density under the HDL 
approximation is smaller than that of the free quark gas. 

So far we have derived an expression of the nontrivial p dependence of the pressure 
density in the framework of the HDL approximation. Now let us give a discussion of the 
features of our approach. The study of EOS of quark matter is a long-standing subject of 
QCD study. Just as was mentioned in the introduction, although lattice QCD has achieved 
some success in treating small p physics, they are far from being complete in the case of 
large chemical potential problem. So one expects that continuum model studies should be 
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complementary to the lattice simulations in the study of QCD thermodynamics at finite /x. 
When the chemical potential is large enough, one may naively expects that perturbative 
calculation is feasible because of asymptotic freedom. However, as was pointed out in Ref. 



26|, the perturbative calculation of partition function of QCD at finite T or /i encounters the 
convergence problem and one has to resort to resummation techniques. From the point of 
view of functional integral approach, the calculation of partition function of QCD amounts 
to the calculation of all possible vacuum bubble diagrams. This is a hard work. In order 
to avoid the problem of calculating the complicated vacuum bubble diagrams, in our work 



we adopt the approach proposed in Ref. 12j, in which the nontrivial fi dependence of the 



pressure density is totally determined by the full quark propagator at finite /i (the quark 
propagator is the simplest Green function of QCD). If the input dressed quark propagator 
at finite fi is exact, the pressure density calculated with the aid of it will also be exact. 
At present it is not possible to determine reliably the dressed quark propagator at finite fi 
from first principles of QCD. So one has to resort to various models inspired by QCD. It 
is generally believed that the HTL/HDL approximations are good approximations of QCD 
when T or /i is large enough. So in the present work we employ the HDL approximation to 
calculate the EOS of quark matter. Compared with other approaches of calculating the EOS 
of quark matter in the literature, our approach is simpler and has a better QCD foundation. 
The behavior of the pressure density under HDL approximation shown in Fig. 3 is 
qualitatively differs from that of the naive perturbation theory. For example, in Ref. |27| 
Fraga, Pisarski and Schaffner-Bielich (EPS) proposed the following cold, perturbative EOS 
of QCD 



!^U-2(^ 



47r^ ^ IT 



a.s /. , 2 ,, \ , A /as 



G^N,^."^^{n-lN,)^.'l{'^y , (31) 



TT V 3 '^ //_ 

where G = Gq - 0.536A^/ + Nf In Nj, Gq = 10.374 ± 0.13 and A is the renormalization 
subtraction point. The scale dependence of the strong coupling constant ^^(A) is taken as 



«s(A) = -— 






where u = ln(AVA|^), f3o = 11 - 2Nf/3, f3i = 51 - 19Nf/3, and (32 = 2857 - 5033iV//9 + 



325N'j/27. For Nj = 3, A^g^ = 365 MeV. The only freedom in the model of Ref. |27| is the 
choice of the ratio A/yU, which is taken to be 2 in that reference. The comparison between 
the HDL result and EPS result is shown in Fig. 3. From Fig. 3 it can be seen that when n is 



smaller than about 1 GeV, the HDL pressure density is smaller than the perturbative result 
to the Us order, but is larger than that to the a1 order. When /x is larger than about 1 GeV, 
the HDL pressure density is smaller than the FPS one, and an important result is that when 
/i tends to infinity, the HDL pressure density tends to the free quark gas result much more 
slowly than does the FPS one. As a comparison, the EOS obtained under Dyson-Schwinger 
equations (DSEs) approach [121] with different parameters are also shown in Fig. 3. 
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FIG. 3: The different EOS 



It should be noted that in drawing the EOSs shown in Fig. 3, we have neglected the 
term 'P(/i)|^=o. This does not mean that the term P(/i)|^=o is unimportant. In fact, it is 
an important quantity. It enters the energy density, which is relevant for integrating the 
Tolman-Oppenheimer-Volkoff equations. At present it is not possible to calculate reliably 
V{^)\^=o from first principles of QCD. One can only calculate it using various models inspired 



by QCD, for example, using the CJT effective action (see, for instance, Ref. [28|). It is 



obvious that the existing calculations of P(/i)|^=o in the literatures are all model dependent. 
In our work, we have not considered the term P(/i)|^=o- The main reason is that this term 
does not affect the relative relation of different EOSs displayed in Fig. 3, since all of the 
curves shown in Fig. 3 do not contain this term. If one tries to add the term P(/i)|^=o in the 
EOS and then make comparison between different EOSs, one has to add it in all different 
EOSs. Therefore the main conclusion of the present paper does not change. 

Now it is time to discuss the range of applicability of the EOS presented in this paper. 
For the low temperature and high chemical potential regime, our knowledge is quite limited. 

10 



Experimentally, it is impossible to achieve such a condition in laboratories on earth. A 
natural laboratory holding such cold, highly compressed matter is a quark star. So the 



only obvious relevance of the presented EOS is for quark stars. It is wel 
astrophysics the study of quark star depends crucially on the assumed EOS 



_cnown that in 
l|, |2|. Therefore 



it is interesting to apply the present EOS to the study of quark star. In our future work 



29| we shall consider the structure of the quark star for the three flavor situation, where we 
shall study the mass-radius relation, mass-central-density relation, and distribution of inner 
mass and chemical potential for this kind of quark star. In addition, we shall compare the 
results calculated from the present EOS with those calculated from other EOSs in Fig. 3. 

To summarize, in this paper we calculate the EOS of QCD at zero temperature and finite 
quark chemical potential under the HDL approximation. It is found that when /i tends to 
infinity the HDL pressure density tends to the free quark gas result much more slowly than 
that of the cold, perturbative EOS of QCD of Fraga, Pisarski and Schaffner-Bielich. We 
also give the expression of QNS at finite T and U-lt is found that when /i — )■ 0, the result 



of QNS is the same as the one obtained in Ref. 25 1. 
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